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We study properties of entangled systems in the (mainly non-relativistic) second quantization 
formalism. This is then applied to interacting and non-interacting bosons and fermions and the 
differences between the two are discussed. We present a general formalism to show how entanglement 
changes with the change of modes of the system. This is illustrated with examples such as the Bose 
condensation and the Unruh effect. It is then shown that a non- interacting collection of fermions 
at zero temperature can be entangled in spin providing that their distances do not exceed the 
inverse Fermi wavenumber. Beyond this distance all bipartite entanglement vanishes, although 
classical correlations still persist. We compute the entanglement of formation as well as the mutual 
information for two spin-correlated electrons as a function of their distance. The analogous non- 
interacting collection of bosons displays no entanglement in the internal degrees of freedom. We 
show how to generalize our analysis of the entanglement in the internal degrees of freedom to an 
arbitrary number of particles. 

I. INTRODUCTION 

Entanglement is a phenomenon that has become of great importance in modern applications of quantum mechanics 
[1]. We usually think of entanglement as existing between different degrees of freedom of two or more particles. There 
has been a number of important advances recently in understanding entanglement in systems containing a small 
number of particles. If we are to use this appropriately, however, we need to extend our analysis to realistic systems 
containing a large number of particles. When it comes to large systems in quantum mechanics, the concept of a 
particle actually fades away and is replaced by the notion of "an excitation of a given mode of the field representing 
the particle". Individual particles actually become indistinguishable (this is, of course, also true for a small number 
of identical particles, but indistinguishability may not play any significant role). In addition to that, the concepts 
of different particle statistics (fermions versus bosons) then also becomes directly relevant and it is important to 
understand its relation to entanglement. The most convenient and appropriate formalism to deal with all the issues 
involving a large number of systems is the second quantization. So we need to understand what entanglement means 
in this setting if we are to be able to harness solid state and condensed matter systems for information processing 
purposes. Another benefit of the second quantization is that it is also the correct formalism for the (relativistic) 
quantum field theory. We can therefore exploit some already known results in this field in the hope of reaching a 
more complete understanding of entanglement. 

In this paper we will see that entanglement is highly dependent on the choice of modes [2] (for an alternative 
treatment see [3]). The choice of modes is mainly dictated by the physics of the given situation, and the general rule 
is that in order for modes to be entangled they need to interact in some way. The crux is that what happens to be 
an interacting hamiltonian for one set of modes may not be so for a different set of modes. Alternatively stated, the 
vacuum state of one set of modes, may not be the vacuum state when the modes are changed. It is this simple fact 
that will be discussed in this paper and shown to be at the root of many different phenomena. More interestingly, 
however, we will see that the modes need not interact in order to produce entanglement - entanglement can be a 
consequence of indistinguishability of quantum particles and the existence of quantum statistics. 

In the next section we review the formalism of second quantization. In section 3 we present a very simple illustrative 
example of two harmonic oscillators and show how entanglement arises through their interaction. The conclusions of 
this section are very general and will then be applied to the Bose condensation as well as the Unruh effect in section 
4. In section 5 we show that interaction is, however, not necessary to generate entanglement. We present a Fermi 
gas of non-interacting particles at zero temperature and show that the spin degree of freedom of individual fermions 
can be entangled even though the fermions themselves are non-interacting. We show that an analogous collection of 
non-interacting bosons, on the other hand, cannot generate any such entanglement in the internal degrees of freedom. 
We finally conclude and discuss some future directions in section 6. 



*e-mail: v.vedral@ic.ac.uk 



1 



II. SECOND QUANTIZATION 



First we briefly review the formalism of the second quantization [4]. We will approach the second quantization 

formalism from the perspective of many particle systems, rather than as a means of achieving rclativistic quantum 
theory. Suppose we have the Schrodinger equation for a single quantum system in one dimension (for simplicity): 

H{x)^{x,t)=ih^^{x,t) (1) 

Note that although this will lead to the non-relativistic second quantization, much of what we say will be true 
relativistically as well. A formal solution to this equation is 

<b{x,t) = J2bn{t)Mx)- (2) 
n 

As usual, by substituting this back into the Schrodinger equation, we obtain its time independent version 

H{x)i>r^{x) = E„Mx) , (3) 

where are the corresponding energies. In order to convert this into a many-particle equation, we apply the 
formalism of second quantization, which means that we effectively have to "upgrade" the wave-function \l> into an 
operator. Formally, we write 

^{x,t) =Y,bn{t)Mx) (4) 

n 

The bji operators are the well-known annihilation operators. The conjugate of the above equation becomes 

4>^X,t)=^bi{t)rn{x) (5) 
n 

where fejj are the creation operators. In this paper we only consider the time independent creation and annihilation 
field operators. Nothing much would change conceptually if the field operators were to be time dependent, although 
mathematically the whole analysis would become much more involved. The Hamiltonian also has to be second 
quantized, and the new Hamiltonian is given by the average of the old (first quantized) Hamiltonian 

H = j dxif''{x,t)H{x)'^{x,t) (6) 

By invoking the orthogonality rules, (V'nIV'm) = ^mn, the second quatized Hamiltonian becomes 

H = Y,Enbik (7) 

n 

and this is the same as a set of independent harmonic oscillators. Now, the fields here are clearly non-interacting and 
therefore we may conclude that the b modes are disentangled. But, we have to be careful. While there is, clearly, no 
entanglement in the b modes, there may be entanglement between some other modes. To illustrate this let us have 
a look at the simplest case of two harmonic oscillators. Note that the index n can also contain the internal degrees 
of freedom (such as spin or polarization), however, not until section 5 will we need to take this complication into 
account. 



III. INTERACTING HARMONIC OSCILLATORS 



Suppose that we have two non-interacting harmonic oscillators so that the total Hamiltonian is (we will omit hats 
from now on): 

H = n{oj-X)b\bi+n{uj + X)b\b2 (8) 

Since oscillators are non-interacting their eigenstates are, therefore, composed of disentangled direct products of the 
eigenstates of the individual harmonic oscillators. For example, the ground state is just a product of the ground states 
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of the oscillator 1 and the oscillator 2, |0)6 = \Q)bi (S) \0)b2J so that it is annihilated by both the annihilation operators 
of individual oscillators 



bi\0)b^O (9) 
b2\0)b = 0. (10) 

However, let us make a change of modes and see what happens to the new eigenstates. Let the new creation and 
annihilation operators be 

(We could have made a more general change, but the one above is sufficient to illustrate the point here). The original 
Hamiltonian now becomes 

H = ?iw{a\ai + 0302) + UX{a\a2 + a^ai) (13) 

and it is transparent that the a modes actually represent two coupled harmonic oscillators. It is, therefore, not 
surprising that the ground state of this system is actually entangled when written in the basis of the original harmonic 
oscillators. It is given by (up to a normalization which is not relevant for our discussion) [5] 

|0)a cx V[-itanh— ]"|n)6, O \n)b, , (14) 

n 

where \n)b. is the state containing n quanta in the oscillator 6j. This state is clearly entangled with respect to 
eigenstates of b s, and it is written in the Schmidt decomposition form making it easy to detect entanglement. The 

amount of (uitanglement is obtained by tracing out one oscillator and calculating the entropy of the resulting mixture 
[1], which is known as the entropy of entanglement. The resulting entanglement is 

21n(coth^) , n4A, 
E= \ ^_ln(l-tanh2— ) . (15) 

To make the result look more familiar to us (and to prepare ourselves for the next section) we introduce the following 
change of variables 

[tanh— ]2=e-tf, (16) 

UJ 

where we think of T as the corresponding (effective) temperature; the entanglement is now given by 

E=^^^-H^-e-'^) (17) 

Ki ekT — I 

which is just the entropy of a thermal state of a harmonic oscillator, i.e. it is the familiar Planck distribution. We see 
that if T — » (i.e. the interaction tends to zero, A — > 0), then E ^ 0, and entanglement disappears as expected. The 
fact that this is a thermal spectrum features strongly in the Unruh effect and the Hawking radiation from a black 
hole [6] as will be seen in the next section. 

Let us now generalize our consideration to an arbitrary field which can be written as superposition of the creation 
and annihilations operators for the various modes 

^{x) = ^(^„(a;)a„ + C(a;)4) ■ (18) 

n 

Here we basically have infinitely many harmonic oscillators to sum over. The vacuum state for the a operators is 
defined as a„|0) = for all n. Now, suppose that we expand the field in terms of a different set of creation and 
annihilation operators, like so 
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n 



(19) 



The new operators can always be expanded as a combination of the old operators as 

&m = ^ Oimnan + l^mnol (20) 

n 

am = Yl a*ranK " Pm-nPl, (21) 
n 

(the new mode functions V' can also be expressed in terms of the old ones [6], but this is not important for our 
discussion). Due to the unitarity requirement, the coefficients of transformation have to posses the following properties: 

Efc (a^fcaJ*fc - /3jfc/3jfc) = h (22) 
X^fc {cxikPjk - Pikajk) = . (23) 

Now, the point of our argument here is that the new vacuum, defined as bm\0) = (for all m), will not in general be 
the same as the old vacuum (i.e. it will not be annihilated by the old annihilation operators, a„). The new vacuum 
will, for example, have the following number of particles 

(6|al«,|0)=^|A,f (24) 
j 

with respect to the ith mode of a (and this means that the vacuum of the b modes contains some particles with 

respect to the a modes). Therefore, since the vacuum |0) is an overall pure state, we can conclude that the amount 
of entanglement between the ith mode and all the other modes taken together is given by the von Neumman entropy 
of entanglement [1] 

Ei = -Y,\Pij\'^^\/3ij\' (25) 
j 

We see that this entanglement is zero only when one of the /3s is a unity and the rest are zero, i.e. the new modes are 
the same as the old ones. The maximum entanglement, on the other hand, is achieved when all the betas are equal. 
Our two modes at the beginning of the section are a very special case of this analysis. 

In summary, there are two main messages in this section: one is that entanglement is highly dependent on the way 
we decide to divide the whole system into different modes; two is that interaction in general gives rise to entanglement. 
There are numerous examples of various kinds of mode transformations in quantum optics (for example, in two-mode 
squeezing) , condensed matter physics as well as relativistic quantum field theory which lead to entanglement generation 
as described above. We discuss a few of them next. 



IV. BOSE CONDENSATES AND THE UNRUH EFFECT 



We first briefly review the second quantization treatment of interacting systems and then apply it to several cases. 
The Hamiltonian in the first quantized form is now given by 

H{x, x') = Ho{x) + V{x, x') (26) 

where Hq{x) is the basic (free) Hamiltonian without interactions and the interaction itself is encoded in the term 
V{x,x'). We already know how to second quantize Hq{x), and so we only need to be able to second quantize the 
interaction Hamiltonian. We will do this now with the intention of obtaining a typical Hamiltonian for the Bose- 
Einstein gas. This is done by performing the following transformation 

V = j dV^\x,t)^\x',t)V{x,x')^{x,t)^{x',t) . (27) 

As before, wc need to switch from the creation and annihilation operators in the position picture to the creation and 
annihilation operators in the momentum picture. This is performed in the same way as we did for the non-interacting 
fields, as in eq. (4). After some manipulations and approximations we can arrive at [7] 
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2 

^ = J2 ^^"'l^'P + S ^(9)«Pl+q42-9«P2«Pl • (28) 
P Pl>P2,9 

So, the first term of the Hamihonian is (as before) a collection of harmonic oscillators and the second term signifies 
their interaction. The interaction term tells us that two particles of momenta pi and p2 collide to produce new particles 
with momenta pi + q and P2 — q respectively. The values of the final momenta are such as to obey the momentum 
conservation law. V{q) is the strength of interaction which is derived from the original potential V{x, x'). The exact 
form of this is unimportant. In fact, it is sometimes customary to assume that the interaction is independent of the 
momentum change q and to write (after omitting the tilde) 

^=ll^'^l'^P+^ «P1+X2-9«P2«P1 • (29) 

P PliP2,9 

This becomes the same as the Hamiltonian for a degenerate, but almost ideal Bose gas [7]. The creation and 
annihilations operators, of course, obey the usual bosonic commutation relations. In order to solve and find the 
eigenvalues of this Hamiltonian we need to convert it into a bunch of non- interacting harmonic oscillators (or modes). 
The transformations that achieve this are called the Bogoliubov transformations and are given by 

Op = Upbp + Vpb^_p (30) 
p = Upb^ + Vpb^p (31) 



a 



where the coefficients u and v (functions of the momentum p) have to be chosen so that the Hamiltonian for b 
modes has no interacting parts (see [8] for an experimental confirmation of the Bogoliubov transformations in a Bose 
condensate). The resulting Hamiltonian, written in terms of the b modes, is given by 

H = Eo + Y,E{p)blbp (32) 

where Eq is, for us an unimportant, constant and E{p) is the energy of new modes (all the details, which are irrelevant 
for the analysis here, can be found in [7]). Here we have a collection of decoupled (i.e. non- interacting) modes as 
in the previous section when we discussed the example of two harmonic oscillators. Again, the ground state is an 
entangled state in terms of the a modes and is given by 

l^o) = n - E{-r>>-^ = t;np = ^. (33) 

In this state only the p and —p modes are entangled and different momenta are not coupled. The amount of 
entanglement is again easy to compute to be (c.f. eq. (17)) 

^=i:((^;:^-Mi-K/»,«) (34) 

This is a sum of entanglements in each individual pair of modes with the opposite momenta. The entanglement is 

highest when Up « Vp. This happens in the low momentum regime, more precisely when p « mc. In the other 
regime, when the momenta are high we have that Up « 1 and Vp ^ and so entanglement drops to zero. 

A very closely related, but even more intriguing, case of mode transformations generating entanglement can be 
found in the domain of relativistic quantum field theory. Next we discuss the accelerated vacuum leading to the 
Unruh effect [6]. The key point in this effect is that the vacuum state for an uniformly moving observer in the 
flat Minkowski space, which is a continuum of decoupled, disentangled, harmonic oscillators, is actually seen by an 
accelerated (Rindler) observer to be entangled. The entanglement is created because the creation operator for a 
particle of momentum p for the accelerated observer needs to be expanded as a sum of the creation operations of the 
uniform observer involving both the p and —p momenta. More precisely, Minkowski quanta (i.e. those detected by an 
inertial observer) are obtained by quantizing the positive frequency solutions of the Klein-Gordon equation; however, 
once we change the time parameter which is what happens for the accelerated Rindler observer, we obtain both the 
set of positive as well as the negative frequency solutions. These, when quantized, lead to completely different set 
of raising and lowering operators, and therefore a new type of quanta. Interestingly, the transformation of the mode 
operators between the stationary and the accelerated observer is actually of the Bogoliubov type discussed already in 
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the case of a Bose condensate. The vacuum state of the accelerated observer is, therefore, an entangled superposition 
of the number states for p and —p of the original uniform field [9] . The interesting observation is that the spectrum of 
the reduced density matrix (when we trace out the negative momenta for example) is Planckian (just like in the case 
of two interacting harmonic oscillators) whose temperature is given by T = a/27rfc, where a is the acceleration of the 
observer and k is the Boltzmann constant. There is, therefore, entanglement between positive and negative modes due 
to acceleration, and this entanglement is proportional to the amount of acceleration. The amount of entanglement, 
more precisely, is the same as in eq. (17) with the appropriate value for the temperature. 

We would just like to close this section by briefly mentioning an effect closely related to the entanglement in the 
accelerated vaciium and the corresponding thermal bath. This effect was discovered by Hawking and it explains the 
radiation from a black hole and the fact that we can associate a temperature and entropy with it [10]. The two effects 
- accelerated observers and the black hole radiation - are, loosely speaking, related because of Einstein's equivalence 
principle, which stipulates that an accelerated observer is the same as an observer in a gravitational field whose 
strength is that same as the value of acceleration. So, someone sitting close to a black hole for example, should also 
see a radiation (just like an accelerated observer does) which has a black body spectrum. The temperature of the 
radiation can be shown to be T = njli^k, where k is the surface gravity. Hawking [10] showed in 1974 that this 
exactly is the case. The entanglement here is generated because the incoming modes which contain only positive 
frequency waves will be converted into the outgoing waves, which will depend also on the negative frequency waves 
due to the metric change during the collapse of matter to form a black hole (there is a pair creation of this type under 
more general metric changes in the expanding universe as first discussed by Parker in [11]. He found that particle 
creation occurs because the positive and negative parts of the field become mixed during the expansion, such that 
the creation and annihilation operators at one time t are a linear (Bogoliubov) combination of those at an earlier 
time.). So, the situation here can mathematically be described by the same Bogoliubov transformation discussed 
before [6,10]. Interestingly, the entropy of a black hole can also be seen to arise due to entanglement between the 
modes outside the black hole and modes inside the black hole [12] (this is a static, rather than a dynamic explanation 
of the entropy of the black hole). There are many issues involved here which have not been resolved, but which are 
beyond the motivation of the current paper, and I direct the interested reader to consult Wald's excellent review on 
this and related subjects [13]. 

The above examples illustrate the entanglement properties in the second quantization and its dependence of the 
choice of modes. A collection of harmonic oscillators that does not interact becomes interacting when we change to 
a different set of modes. This change may be "passive" , like in the case when we diagonalize the Bose condensate 
Hamiltonian, or it could be "active", like when a stationary observer accelerates. But, the question we would really 
like to ask in this paper is whether entanglement could arise without any interaction between different modes. So far, 
we have only considered spinless bosons, and the particle statistics has not been crucial in any respect. We would 
now like to analyze the interplay between the internal and external degrees of freedom in fermions and bosons, which 
are otherwise non-interacting. 



V. NON-INTERACTING FERMIONS AND BOSONS 



Suppose we have a collection of non-interacting electrons, in a metal, for example. The electrons all move with 
various different momenta p and Pauli's exclusion principle tells us that there can be only two electrons with the 
same value of momentum (one with spin, s, pointing "up" and the other with the spin pointing "down"). The state 
of these electrons is well approximated by the overall state (the so called Fermi model of a metal): 

l'Ao)=n^lW|0)- (35) 

The fact that these are Fermions (as opposed to Bosons) is reflected in the anticommutation relations which we have 
to impose between the annihilation (and creation) operators 

\h\{p)M{<l)\+=^st^{v-<l) . (36) 

Obviously there arc no correlations between the momenta of electrons. This means that if we find an electron and 
measure its momentum to be p, then if we measure another electron's momentum its value will be completely 
independent of the first measurement. That is exactly what the product of creation operators in the above state 
signifies. However, there is a great deal of entanglement in the above state, we just need to look for it more carefully. 
Suppose that we measure an electron at the position r and another one at the position r' . We will assume, for the 
sake of simplicity, that we have arbitrarily precise detectors (this is not a point of fundamental importance in the 
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analysis). Suppose further that one electron has spin up, can wc infer the direction of the spin of the other electron? 
The question is, are electrons at different locations correlated in in their spin directions? In order to calculate this, 
we need to compute the density matrix describing the spin state of two electrons. It is (by definition) given by 

Pss';u' = {<Po\4'{r')4{r)^Ar')i^s{r)\M (37) 

where il)\{r) creates a particle of spin t at the location r. This is the same as computing the elements \st){s't'\ of the 
two-electron density matrix. Assuming (as we usually do) that 

^s{r) = J ^e''"'bs{p) (38) 
we obtain the following form of the density matrix 

Pss ;tt'-J (27^)3 (27r)3(27r)3(27r)3'' 

X {MbUp)bl{q)bs'{q')bs{p')\<l>o) ■ (39) 

Now we evaluate the last term to be 

{Mbl{p)bl{q)bs'{q'%{p'Mo) = {27TfS,tSs't'S{p-q)6{p' - q') 

- {27rfSs'tSst'S{p' - q)S{p - q') . (40) 

This form is intuitively clear as we have both the direct term, i.e. the term where the particle with momentum p 
changes to q and p' to g', and the exchange term, such that the particle with momentum p changes to q' and p' to 
q. Note that because of the conservation of momenta these have to be the same (there is no interaction to allow the 
change of momenta or conversion of one particle into two for example). This now enables us to compute the density 
matrix, by observing that the density of electrons is given by n = J^" dp^ / (27r)^ (where pp is the Fermi momentum), 

Pss';tt' = n^i^stSs't' - 5s't5st'f{r - r')) , (41) 

where 



-r')= / 
Jo 



f(r-r')- I ^e-ip{r-r') ^ 3ji(fcf|r-/|) 

■'^ ^ ' (27r)3 kF\r-r'\ ^ ' 



is the "exchange interaction" term, and ji{kF\r — r'\) is the spherical Bessel function with kp = Pr/fi- This form of 
p contains all 16 elements needed for the spin density matrix. This 4 by 4 density matrix is (up to a normalization) 

given by 



Pl2 



/ 1-/2 \ 

1-/^0 

0-/2 1 

V 1 -/2/ 



where the subscript 12 signifies that there are two electrons, and the basis is the usual | It), | Ti)) I it)) I ii)- We first 
have to decide if this matrix is entangled. We perform the partial transposition and compute the resulting eigenvalues 
[14]. They are 1-/2 (doubly degenerate), 1 and 1 — 2/2. The first two eigenvalues are positive as < / < 1, but 
the last one can be negative. The condition is that /2 > 1/2, and in this region we have entanglement. Therefore the 
density matrix is entangled for the region of distances < |r — r'| < Vg, where is the value such that jf{kpr,,) = 1/2 
(the value is about ~ nh/SpF as expressed in terms of Fermi's momentum). 

We can compute the entanglement of formation [15] of the spins for a general state pi2 (as a function of /). It is 
given by 



E{pi2) = ^ In 



■ 2— f 2 , 

In ' ' ^ (43) 
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Entanglement is maximal when / = 1 and it is equal to In 2. The resulting state is the spin singlet state | "fj,) — | ||). 
Otherwise, as the distance decreases and / drops, the more of the triplet state becomes mixed in with the singlet 
and entanglement drops as a consequence. Beyond the distances of about nh/Spp, all entanglement disappears. 
A reasonable estimate for the Fermi momentum in a metal is 10~^'^ kglO'^ms~^ , which gives a distance of about 
Ar > 10~^ meters beyond which no two electrons are entangled in spin. Below this distance, every two electrons in 
a conductor are entangled. If we have a much bigger momentum, such as in a neutron star, then the entanglement 
distance becomes very small. For a typical neutron star, this distance would be about 10~^^ meters, which tells us 
that there is no entanglement there between any two neutrons in a neutron star. Therefore, the higher the Fermi 
momentum the smaller the "entanglement distance" between the electrons. Note that it is a separate issue whether 
this entanglement present in the spin degrees of freedom can actually be used to process information, such as in the 
teleportation protocol. In order to teleport, for example, the electrons need to be separated sufficiently well after 
being detected. How feasible this is in practice will be a subject of a separate investigation. 

The disappearance of entanglement is intuitively expected from the following simple argument. Quantum statistics 
becomes important only when the spatial wavefunctions of individual fermions have a substantial overlap. The 
wavefunction spread of fermions is roughly given by the wavelength which can be estimated to be h/pp- Therefore 
beyond this distance we expect particle statistics not to play such a strong role. Let us look at this in the first 
quantization to make the case even clearer. Suppose that the total (un-normalized) state of two electrons is: 

|01)|Vi,V2) (44) 

where the {ipi, '(p2) are the spatial wavefunctions of the two electrons, and 0, 1 are the values of spin (the first position 
always belongs to the first electron and the second to the second electron). Suppose that ('!/'i|V'2) = e- Now this is not 
an allowed state if the particles are truly indistinguishable and in the case of fermions it needs to be antisymmetrized. 
The corresponding allowed state is 

|01)|Vi,V'2)-|10)|V2,^i) (45) 

Now, this state can be entangled in spin. To check this, we trace out the spatial degree of freedom to obtain the 
following density matrix (again unnormalised) 

|01)(01| + |10)(10|-e2(|01)(10| + |10)(01|) (46) 

It can be checked, again using the partial transposition criterion, that this state is entangled, unless e = 0, which 
means that the two electrons have no spatial overlap. Therefore we see that the quantum statistics plays a crucial 
role for entanglement. 

Classical correlations, on the other hand, are still present even beyond this distance. This is true in the above first 
quantized model, as well as in the model of many non-interacting fermions. We can quantify them using the mutual 

information, I{pi2) = S{pi) + S{p2) — S{pi2), in the state in eq. (41) [1]. Here pi and p2 are the reduced density 
matrices derived from pi2. Computing this for our state of two electrons we obtain 

/(p,,) = 21n2 + 2^1n^ 

1 1 1-2/^ 1-2/^ 

^ 4 - 2/2 4 - 2/2 + 4 - 2/2 4 - 2/2 ^ ' 

We see that the classical correlations only vanish when / = as we then have an equal mixture of all four states 
I T) T)) I T) i)) I i) T)) I i) i)- Otherwise, spins are always (classically) correlated. 

The entanglement in spin when space is traced out is of a similar (but somewhat dual) nature to the spin-to-space 
entanglement transfer in the double beam-splitter set up in [16]. The entanglement, which in this experiment by Omar 
el. al. initially exists in the spin degrees of freedom, is, due to the quantum statistical nature of particle interaction, 
transferred to space. This also is a purely quantum statistical effect as there is no other interaction between the 
internal and external degrees of freedom (although there is an interaction between the spatial degrees of freedom of 
different system involved). 

Note that the presented method here allows us to extract more complicated information such as correlations between 

three and more electrons. These would be increasingly more difficult to treat mathematically, but the calculations 
can be performed in principle. We could, for example, check if genuine multi-electron correlations exist and what are 
the relevant distances to allow such correlations. Let us briefiy summarize the calculation: we have three electrons 
at positions r, r'r" and they have some values of spin. We would like to extract the spin density matrix, which is 
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now at 8 by 8 square matrix. The form of the matrix can actually be deduced without much effort from the simple 
"pictorial" representation. We have the following terms 



f'^-y 1 1' t" [ 1 1" t' ] [t'tt" [ t" t' t r' 



s 



s s' s" 

t' t" t 



s ,s ,s 

t" t t' 



There are six terms as there are six possible arrangements; every spin in the density matrix can initially be paired 
with any of the three spins. The next spin has only two possibilities and the last one has only one remaining choice. 
The choice of signs should also be clear: the first term is the direct term and therefore positive. The next three 
terms involve one double exchange and arc therefore negative, and the last two terms involve triple exchanges and 
are therefore positive. The density is, therefore, given by 



p{s, s' , s"; t, t' , t" 



= (Mi^l, {r")4 ir')4 {r)7p, (r')V«' (OV-^" {r"Mo) 

+ v?{5si5s'i,'5s"t" 

- f2f3Sst"Ss't'5s"t 

+ fif2h5st'5s't"5s"t 
+ hf2h5st"5s't5s"t') 



(48) 



3/2 ^ 


-2/3 
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There are three different / fimctions now as we have three arguments |r — r'|, |r — r'| and \r' — r"\. This contains 
information about all 64 elements. The actual matrix is (for simplicity in notation, I will assume that all the / 
functions are equal) 

/I 



Pl23 



V 

In order to confirm that we obtain the same 2 electron density matrix, we have to trace one electron out, and this 
we have to do by tracing over both the spin and space. If this is performed, we can easily confirm that the resulting 
density matrix is the one obtained previously for two electrons. 

Now, we would like to know if this state contains three party entanglement in spin. We already know that every 
two parties can be entangled in the appropriate distance region, for example, when they are very close to each other. 
What happens when r = r' = r"l Then, fi = l and we see that the resulting state is actually not a physical state 
(all the elements go to zero indicating that this is not a physical state - or that the probability for this to happen 
vanishes). This is physically clear, as when we have three spins, at least two out of three have to be in the same state 
and this is not possible as they are all in the same physical location. At the other extreme, when all three electrons 
are very far away, then the density matrix again becomes maximally mixed. There is no entanglement in this case of 
any kind, nor do we have classical correlations. We also have a number of other possibilities, such as, for example, 
that two electrons are close to each other, but the third one is further away from them. We could then imagine that 
all three are in an entangled state as well as having entanglement between the two close by electrons. Let us take as 
an example the state when all electrons are very close to each other. Then, to a good approximation, the state is an 
equal mixture of the four state |000), |111), |001) + |010) + |100) and |011) + |110) + |101). Now, this state clearly 
contains three party entanglement. This is because every electron is entangled to the other two electrons, indicating 
that the state cannot be separable in any way. Following this logic therefore we have a three party entanglement 
providing p > 1/2. Can we have a genuine three part entanglement beyond this point even though there is no two 
party entanglement? This would be something like a GHZ state (where there is no two party entanglement, but all 
three parties are entangled) and this remains an open question. Generalizing our line of reasoning, however, we can 
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state that if a bunch of n electrons is found within the distance of roughly the inverse Fermi wavenumber, then the 
corresponding total state of their spins would be an entangled state on all levels (two, three-party and so on up to 
n-party entanglement). The above treatment also allows us to compute various mutual informations between different 
parties involved. 

Interestingly, an equivalent analysis can be performed for a bunch of non-interacting bosons (with polarization as 
the internal degree, for example), but we will not see any correlations in polarization. This is because there is no 
restriction on how many bosons can be in a state with a given momentum and polarization. The most general state 
of bosons is now given by 

\<l>}^U(blip)r^'''^\0} (49) 

sp 

Where n{p, s) is the number of bosons of momentum p and polarization s. The operators obey the following commu- 
tation relations 

[bl{p),bt{q)]-=6,t5{p-q) (50) 

It is now intuitive that if we measure a photon of one polarization at r, there is no reason why this should be correlated 
to the polarization of another photon at r'. We can compute the same density matrix as in the case of fermions, 
but the quantity {(l>o\bl,{p)bl{q)bs'{q')bs{p')\(po) will now be independent of the labels s,s',t,t'. This means that the 
resulting 4 by 4 density matrix of polarization will have all the 16 elements equal to each other. This represents 
a disentangled pure state {\H) + \ V)) (g) {\H) + \V)), where H and V are the horizontal and vertical polarizations 
respectively (the polarization plays the same role here as the spin did in the case of fermions). And, so, here we 
see a fundamental difference between fermions and bosons with respect to the interplay between the entanglement in 
internal and spatial degrees of freedom. Therefore, there is no correlation between the relative distance between the 
photons and their polarization. Bosons, in other words, do not obey the Pauli exclusion principle. 

Note that in this section we have analyzed the presence of entanglement in the internal degrees of freedom and have 
not talked about the spatial degrees of freedom. The analysis of long range correlations in fermions and bosons can 
be found in [17], although the author there does not distinguish between classical and quantum correlations (and the 
internal degrees of freedom play no special part). Nevertheless, similar methods could possibly be used in the future 
to perform a more detailed analysis of more complicated systems. Note also that all our discussions here could be 
generalized to overall mixed state (here the total states we start with are pure states) as in the statistical mechanics 
of many body systems [18]. The same analysis can be performed as the one presented above to calculate the amount 
of entanglement in the internal degrees of freedom. It is to be expected, however, that this additional mixedness can 
only reduce the existing correlations in an overall pure state computed here. Otherwise, the results would not be 
fundamentally different to the ones presented here. 

VI. CONCLUDING REMARKS 

We have seen that the second quantization formalism adds some interesting effects to our understanding of entan- 
glement. First and foremost, this is because we have to decide what our subsystems arc before wc can speak aboiit 
the amount of entanglement between them. This is not an issue in the first quantized version of quantum mechanics 
that is used when we deal with a small munber of systems, such as 2 or 3 (by definition) distinguishable) qubits. An 
interesting observation is that the change of modes also affects the amount of entanglement present. A change of 
modes is, of course, a "global change" relating the old modes to the new ones and this is why entanglement can be 
affected in the first place - "local changes" on their own would not be able to generate entanglement (We have to be 
a bit more careful here as we have to restrict ourselves to speaking about modes and excitations of the modes. We 
cannot really speak about local operations on "particles" , since they are identical and whatever operation we perform 
on one of them, we, in fact, have to perform on all of them). The change of modes produces an interaction between 
the new modes which is what gives rise to entanglement. In section 3, we presented the most general formalism for 
mode changes and showed how the amount of entanglement varies with the change of modes. However, we have also 
been able to show that we do not really need an interaction to produce entanglement. We can instead rely on the 
effects of quantum particle statistics to produce entanglement. This effect was illustrated by looking at the Fermi sea 
model of a conductor, where we had a collection of non-interacting electrons (neither in spin nor in space) occupying 
states up to the highest Fermi level (assumed arbitrary). Here we found that two electrons are entangled in the spin 
degrees of freedom providing that their distance is (roughly) below the inverse Fermi momentum (measured in units of 
fi= 1). Our method can be used to compute the density matrix for any number of subsystems. We have also pointed 
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out a crucial difference between non-interacting fermions and non-interacting bosons, in that non-interacting fermions 
can still be found to be entangled in the internal degrees of freedom, while bosons cannot. It would be interesting to 
investigate the further possibility of using these "non- interacting" correlations to perform information processing [19]. 
It would also be interesting to extend our discussion to more complicated systems composed of fermions and bosons. 
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